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Abstract 

We establish an hitherto hidden connection between zero modes and instantons in the context 
of the Kraichnan model for passive scalar turbulent advection, that relies on the hypothesis that 
the production of strong gradients of the scalar is associated with Bose-like condensation of La- 
grangian particles. It opens the way to the computation of scaling exponents of the iV-th order 
structure functions of the scalar by techniques borrowed from many-body theory. To lowest order of 
approximation, scaling exponents are found to increase asymptotically as log N in two dimensions. 

PACS numbers: 47.27.eb, 02.50.-r 
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There has been in the last two decades great progress in the quantitative understanding of 
anomalous scaling in fluid turbulence. In the framework of the Kraichnan model for passive 
scalar advection, the emergence of anomalous scaling has been traced to the existence of 
statistical integrals of motion showing up in the evolution of Lagrangian fluid particles 
(see (l[ for a review). For any finite number of particles, the conserved quantities are 
functions of the interparticle separations that are statistically preserved as the particles are 
transported by the random flow and scale as a power law with respect to the mean radius 
of the cloud of particles. The scaling exponent of these so-called zero modes (for reasons to 
get obvious below) depend in a nonlinear way on the number N of particles and it turns 
out that the scaling behaviour of the N-th order structure functions of the advected scalar 
^iv( r ) = + r ) — @i x )) N ) (f° r N even) is dominated by the irreducible zero mode of 

corresponding order and lowest positive scaling dimension. Although some specific features 
of the Kraichnan model are needed to establish on a firm mathematical ground the existence 
of zero modes, this mechanism is believed to be robust and relevant for transport by generic 
turbulent flows. 

On the other hand, for large values of N, one generally expects iV-point structure func- 
tions to be controlled by rare events which can be captured by the instanton formalism first 



introduced in the field of particle physics (see, e.g., Coleman {2]) and adapted later on for 

n 

turbulence by Falkovich, Migdal and co-workers |3|, at about the same time as the zero 
modes breakthrough took place. In this second approach, one looks for configurations of 
the advecting velocity field (the random noise in the Kraichnan model) of optimal statisti- 
cal weight leading to a prescribed value of the scalar increment at short scales. When the 
dynamics in the inertial range is scale invariant, one may focus on self-similar instantons as 
shown in the framework of shell models of turbulence % |5|- Those objects yield the best 
picture of the singular scaling fluctuations S r 9(x) = 9(x + r) —9{x) ~ r h ( x \ which are at the 
basis of the phenomenological multifractal description of turbulence introduced by Parisi 
and Frisch sometime ago |6(. In the multifractal model, the probability of a h— fluctuation 
occurring is given by P r (h) ~ r s ( h \ where the function s(h) can be interpreted as the extinc- 
tion rate of the singular fluctuation as it cascades towards small scales. Structure functions 
of any order may be computed by averaging over the probability P r (h) the appropriate power 
of the fluctuating field. It is then easily recognized that 7jv(r) oc J dhr ( - s ^ +NhS) is a scaling 
function of r in the inertial range, of the form T/v(r) ~ r^ N with (n = min^,[s(/i) + Nh]. 



Scaling exponents (n are therefore linked by a Legendre transform to the function s(h), 
whose instanton formalism is able to produce a first estimate (dressing of the instanton by 
fluctuations has to be considered in order to improve the result). 

So far, the connection between these two descriptions of the origin of intermittency in 
stochastic turbulent systems has remained elusive. Computations of zero modes and their 
scaling exponents in the framework of the Kraichnan model were mostly done using pertur- 
bative methods around limiting values of its parameters for which anomalous scaling disap- 
pears (high space dimensionality, small or large roughness degree of the advecting velocity 
field). Those methods fail to capture the scaling exponents in the large iV nonperturbative 
domain (reached, roughly speaking, when iV gets larger than the space dimensionality d). 
Lagrangian numerical methods as introduced in {t| offer an alternative strategy for comput- 
ing zero modes. But in practice, due to the rapid increase of the time of integration of the 
equations of motion needed to reach the required accuracy as the number N of Lagrangian 
particles increases, they are bound to explore moderate values of statistical orders. 

Alternatively, there has been some attempts at developing the instanton formalism for 



the Kraichnan model 



9j. The most convincing one was carried out by Balkovsky and 
Lebedev [9j. In order to get rid of soft modes (a common source of trouble in path integral 
formulations and their treatment by saddle-point approximations), physically associated 
with slow modulations of the strength and the orientation of the local strain matrix, these 
authors wrote down an effective theory concentrating on the evolution of interparticle dis- 
tances, which could be then solved only by considering that the space dimensionality d is 
large. The sequence of bold simplifying assumptions introduced in 9] leaves little hope for 
understanding along this way the link between zero modes and instantons in the large N 
limit. 

It is the main purpose of this Letter to bridge the gap between these two cornerstones of 
the modern explanation of the emergence of anomalous scaling in turbulent advection. In 
order to do so, we shall rely on the crucial observation that the production of large scalar 
differences at microscopic scales requires correlated motion of Lagrangian particles coming 
from different regions of space. In a quantum analogy, this will be achieved provided parti- 
cles, seen as bosons, condense in the same state. Standard methods developed in condensed 
matter physics for Bose condensation will allow us to reduce in the large N limit the N- 
body problem one is facing in the zero mode approach to an effective one-particle problem 



involving both diffusion and advection in the mean velocity field created by the remaining 
N — 1 particles. Estimates for scaling exponents of zero modes can then be inferred from the 
ground state energy of the Bose condensate. Furthermore, it turns out that the equations 
fixing its wave-function can be mapped exactly onto the ones coming out from an instanton 
analysis applied to the original stochastic equations governing the evolution of the full scalar 
field 8(x,t) in <i-dimensional space. A numerical study of these equations, performed in the 
present work at d — 2, leads to the prediction of scalin g ex ponents (n growing like logiV 
at large N. This is at odds with the common belief j^, 10] that the Cn's should saturate 
at large N, due to the statistical preeminence of fronts carrying finite discontinuities of the 
scalar like in Burgers equation. We do not know whether the logarithmic growth of scaling 
exponents found in our approach is a robust feature which will survive upon the inclusion 
of fluctuations around the Hartree-Fock description of the Bose condensate presented in 
this Letter. We note however the absence of a definitive theoretical argument in favor of 
saturation of the scaling exponents and leave the resolution of this discrepancy open for 
further work. 



In the inertial range (where forcing and dissipation can be neglected), the equation of 
motion of the Kraichnan model 

dt9(r,t)+v(r,t).Ve(r,t)=0, (1) 

just describes passive advection of the scalar quantity 6(r,t) in a random Gaussian incom- 
pressible velocity field v(r, t) whose two-point correlations behave like 

(^(r,t)^-(r',0> = 25(t - i')(A)<% - dy(r - r')), 

with dij{v) = ((d — 1 + ~ C^")- The exponent £, which fixes the way velocity 
differences scale at short distances and ranges between and 2, together with the space 
dimensionality d, are the two physically important parameters of the Kraichnan model. In 
eq. 02J, spatial indices i and j run between 1 and d, and the incompressibility of the velocity 
field is warranted by the following property of the matrix dij(r) : VjG^r) = 0, Vj. 

Let us now consider iV fluid particles, close to each other at initial time t = 0, and 
transported at later times in various realizations of the Kraichnan velocity field. We define 
Vn(ti, r 2 , . . . r N , t) as the PDF of their positions r 1 ,r 2 , . . .r N at time t (to keep notations 
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simple, we skip the reference to the initial positions of Lagrangian particles in the arguments 
of Vn- We shall also use in the following r as a short-hand notation for (rj., r-2, ■ • • , rjv))- in 
the translation-invariant sector, Vn evolves as 

dV N 



dt 



d ij( r nm)V ri y rL I V N = M N V N , (2) 
n<m / 

where particles are labelled by integers n or m and vectors r nm = r m — r n denote their 
relative positions. 

The operator Ai n is of dimension £ — 2 with respect to length rescaling. As a consequence, 
there exist non trivial solutions to the equation MnJ = (or zero modes), that are scaling 
functions of positive dimension ( (i.e., such that f(Xr) = A^/(r)). Since M. n is a self- 
adjoint operator, one easily deduces that the statistical average of zero modes (defined in 
the translation-invariant sector as (/)(£) = / f (rlTVfr, t)d Vi • • -dfir^-i is conserved by the 
dynamics. It was shown (see again the review paper [l| for more details) that zero modes 
are formally present in the A-point equal-time correlation function of the scalar field and 
that the exponent (n can be identified with the smallest possible value of scaling dimensions 
of irreducible zero modes of A4n, «.e., those modes which do not belong to the kernel of 
sub-operators M. v with p < N. 

In order to make the computation of (n more amenable to techniques of many body the- 
ory, we first relax the constraint of translation invariance by adding a supplementary particle 
at position riv+i and write Pat as Pa? (r, t) = JVn( v i — r iV+i> • • • , r N ~ r N+i, t) d d rN+i- The 
new function Vn, as the PDF of relative positions in a cloud of iV + 1 particles, obeys the 
equation $Ezl = M. nVn, where 

~ 1 N 

M-n = — {dij(r N +in) + dij(r N +im) — dij(r nm )} x V r ^ (3) 

n,m=l 

is nothing but A^at+i expressed in the referential frame of the (N + l)-th particle (from now 
on, we shall put without loss of generality r N +i = and note r N +i n = r n ~ r N+i as r n ). 
It follows that zero modes of Mn are also zero modes of A4n- Another advantage of this 
new formulation of the problem lies in the fact that interactions or correlations between the 
particles in the original cloud and the central one will be taken into account exactly, even 
in the approximation schemes we shall be obliged to introduce later on to make progress. 
This is close in spirit to the Bethe-Peierls method in statistical mechanics, known to give a 
better description of local ordering in condensed phases than basic mean field theory. 



We then take advantage of the scaling properties of the operator M. n to switch towards a 
representation incorporating in a natural way the average expansion of length scales implied 
by the dynamics (j2J). We rewrite Vn as TVfct) = t ~&n(l , lni) with 7 = 2 — £. After 
defining new space and time variables p n = r n t 1 and r = hit, the evolution of the PDF 
takes the final form : 

d$ N f ~ A N AcH 



c/r ^ 7 7 

where A^r = Pn^p; * s the generator of scale transforms in the Ad-dimensional space 
of configurations and M.n is now expressed in terms of the rescaled position variables p n . 
Unlike M. n, the operator C N is not self-adjoint. Since 1 Cn = M-n — we deduce that zero 
modes are left eigenvectors of Cn (of eigenvalue — ~ if £ denotes their scaling dimension). A 
more detailed analysis reveals that zero modes lie at the top of a "tower" of left eigenvectors 
of Cn of eigenvalue — - — k where k is an integer (those eigenvectors are obtained as linear 



111]). Corresponding right eigenvectors of 



combinations of the slow modes discovered in 
Cn are of finite norm. It should be noted that zero modes are also right eigenvectors of 
Cn, and as such give rise to a positive part in the spectrum of this operator consisting of 
eigenvalues of the form - + — + k (with k an integer). However this part of the spectrum 
cannot take part in the time evolution of the PDF since it is built from non normalizable 
right eigenvectors, and will be of no concern in the following discussion. 

We are now in a good position for catching an estimate for Cn by variational approach. 
We can see indeed (n as the minimum of the functional (^/\ — ^Cn |$) for any pair of left 
and right irreducible states (\&| and |<&) of unit overlap (^|$). Irreducibility of the right 
state |$) (which is enough to project the whole variational procedure on the desired Hilbert 
space) is easily enforced by the condition J <&(p) d d p n = 0, Wn < N. The latter is indeed just 
a way of stating that the trial state |<&) does not belong to the dual space of the kernels of 
sub-operators A4 P with p < N. We assume Bose condensation of particles in a dumbbell-like 
geometry, and restrict our attention to variational states of the form 

N N 

$(p)=n^(Pn), ^(p)=n^(Pn), (4) 

n=l ra=l 

where both orbitals <pq and ipo are °dd with respect to inversion of coordinates along the 
dumbbell axis and invariant with respect to rotations around this axis. Minimizing (^/\ — 



7>Civ|$) under the constraint (V'qIvq) = 1 leads to the following conditions for y? and ipo : 

- l N (p = fi N (p , -Hn^q = fx N ip , (5) 
where //jy is a Lagrangian multiplier and the one-particle operator l N reads 

l N = d ij {p)V i V j + \Viip) + ^]V, + ~ (6) 

7 7 

with 

Vi(p) = (N-l)j d tJ (p - P ')VMp') MP') d d p'. (7) 

Particles are seen to condense in the state of lowest "energy" fi^ of an effective Hamiltonian 
involving both diffusion and advection in an incompressible velocity field V(p) defined by 
eq. (J7J), which expresses on an averaged way the interactions with other particles and adds 
to the radial velocity component issuing from the continuous rescaling of lengths. Note 
that the presence of a diffusion term in the effective Hamiltonian is essential to ensure 
the existence of non trivial and physically meaningful solutions to eq. (jSJ). We obtain the 
following upper bound for ( N : 

4 0) = 7 - N{N ~ 1} J Mp)Mp') dij(p - p')V m {p)V m {p') d d pd d p'} . (8) 

We get another interesting relation from our theory by treating N as a continuous variable. 
Differentiation of Pn = — (^ol^vl^o) with respect to iV leads to the result ^§- = — jv^T^v" 
where we define Sn=— / ^o(p)' ? /'o(p / ) d^p — p')V 'j(^ (p)V,v9 (p / ) d d pd d p' '. Considering 

Sn as a function of \i rather than N (via /j,n), we can rewrite in the large iV limit the 
expression obtained before for as 

C ] = 1 {N f i N + S(fi N )} with S'(fi N ) = -N. (9) 

In other words, up to the pre-factor 7 converting time scales into spatial ones, the scaling 
exponent (n is nothing but the Legendre transform of the function «S(//). We have therefore 
recovered, at this level of approximation, the phenomenological content of the multifractal 
model of Parisi and Frisch. 



Let us now show how the preceding results also arise from the instanton approach. In 
order to do so, we shall follow the formulation of instanton theory introduced in [4|, |5[ , though 
in the restricted framework of shell models. We go back to eq. ([!]) and first eliminate global 



sweeping effects by adopting a quasi-Lagrangian description. This amounts to defining all 
the fields in a frame whose origin moves with the fluid and transforms eq. (JT|) into 

d t 9(r,t) + [v(r,t) - v(0,£)].W(r,t) = 0, (10) 

upon appropriate redefinition of coordinates and fields. We then go from the Stratonovich 
prescription underlying the above stochastic equation to the Ito convention. According to 
standard rules of stochastic calculus [12J, eq. (fTUj) becomes 

d t + [v(r, t) - v(0, t)\.V9 - dij^ViVjO = 0. (11) 

The Ito form of the equation has the merit of making explicit the mixing of the scalar due 
to short-range fluctuations of the velocity-field, in the form of an effective eddy diffusivity 
varying with the relative length scale in the local frame. Due to this eddy diffusivity, the 
formation of fronts in the scalar field cannot be anymore arbitrarily slowed down as formally 
allowed by eq. ffTOj) . Furthermore, the evolution of the scalar field in this representation 
becomes a true Markov process, which facilitates the writing of its PDF as a path integral 
over various realizations of the noise. Indeed, following for instance [3j, it is rather easy to 
establish that the probability of reaching a prescribed scalar configuration 9f(r) at time tf, 
being given an initial one 9i n (r) at time ti n , can be expressed in terms of the Martin-Siggia- 
Rose path integral J T>9Tp exp —S, with the action 

3 = dt^J d d rip(r, t)\d t 6 - d# (r)V \V 'jO] - ~ J d d r rfVa^r, t)D i:j (r, r') aj (r', t) j , 

(12) 

where we defined -Dy(r, r') = c%(r) + dy(r') — ^(r — r') and a(r, t) = —ip(r,t)'V8(r,t). 
In eq. ( fl2l) . p(r, t) is an auxiliary field conjugated to the physical field 9(r, t), enforcing the 
constraint that the equation of motion (jlip be satisfied at any time along the trajectories 
contributing to the PDF, and the summation on all possible realizations of the stochastic 
velocity field at intermediate times has already been performed. It will be convenient to set 
p(r,t) = —irj(r,t) and to consider 7] as real in the following. 

Extremization of the action with respect to the configurations of both fields 9(r) and 
7^(r) between times U n and tf, leads to the following set of coupled dual equations defining 
extremal (or instantonic) trajectories : 

d t 6 + Viiv, t).Vi9 - dij^ViVjO = 0, (13) 
-d tV - Vi(r, t).Vi7] - d^ViVjr] = 0, (14) 
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where the velocity field V(r, t) of components Vi(r, t) = J d d r'Dij(r, r')a,j(r', t) describes self- 
consistently the flow in the moving frame of the instanton. We look for particular solutions of 
those equations describing the formation of fronts of the scalar 9 and self-similar blowing-up 
of its gradient in finite time. They can be parameterized as 

9(r,t) = (t*-t) a 9(p), r ] (r,t) = (U-t)- a -^(p), (15) 

where we defined as the (arbitrary) critical time of the blowing-up, the scaling position 
variable p as p = r(t* — t)~~ , and introduced a scaling exponent a which is bound to be 
positive due to the conservation of the scalar. The two scaling functions rj(p) and 8(p) are 
assumed to be odd with respect to inversion along the compression axis of the flow. 

By plugging the Ansatz ( fl5l) in eqs. ( |T3l) and f fl^j) . it is straightforward to check that 9, 
rf and a obey the same equations ([5]), (Q and ([7]) as ip , (p and fijy, provided the overlap 
C = (?} | 9) of the scaling functions (which is conserved along extremal trajectories) is iden- 
tified with N — 1. For a given C, the dynamically selected instanton is the one of smallest 
scaling exponent a, and its action 5* takes the expression —S(a) log *** } (or —75(a) log ~ 
if L and r denote typical length scales of the structure at initial and final times), where 
S(a) turns out to be the same function as «S(/x) introduced in eq. (EJ). We conclude that 
there is indeed a complete merging of zero modes and instanton pictures in the large N limit. 

We solved numerically eqs. ([MZ]) for d = 2 and 7 = 4/3. In order to determine po, 
■00 and fijy for a given value of N, we integrate forward in time the couple of equations 
¥ = h\<p) ~ i f^%> and d -f = (H N iP\ - until a fixed point is reached. The 

equations were discretized within a rectangular box \x\ < L x , \y\ < L y , with L y = 20 and 
L x ranging between 20 and 40 according to the value of N (we assume that the particles 
collapse onto the x-axis, so that the increase of typical lengths with iV is more pronounced 
along that axis than along transverse directions), and a uniform mesh size Ax = Ay = 0.2. 
Some care has to be exercised in the discretizing of fields and differential operators, in order 
to get a sensible translation of the condition of incompressible flow on a lattice. The initial 
conditions for ip(x,y) and ip(x,y) are arbitrary, except for being odd (resp. even) functions 
of x (resp. y). As to boundary conditions, we impose the vanishing of if and the normal 
derivative of ip on the perimeter of the box. We followed the evolution of solutions up to very 
large values of iV (of the order of 10 6 ). Figured] depicts their shape along the positive x-axis 
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FIG. 1. Plot of the one-particle orbitals ifQ and tpo along the instanton symmetry axis for N = 
4, 8 and 16. Curves go from top to bottom and left to right as iV increases. We normalized tpQ by 
setting its slope at the origin to the constant value 0.04. 

for N = 4, 8 and 16. It is seen that bending of the orbital i[)q (the scalar field configuration 
in the instanton picture) goes together with an enlargement of the range of the localized 
dual orbital ifQ. This remains true at higher values of N. 

The one-particle energy fi N was found to decrease asymptotically as l/N, so that the 
action Sn, from eqs. (jSJ) and (jUJ), behaves as logiV and provides the dominant contribution 
to the scaling exponent (n at large values of N (see fig. [2]). The linear extension of the 
dumbbell formed by the cloud of particles, as measured by the distance between the two 
maxima of \(po{x,y)\ along the x-axis, turns out to scale also like logiV, while its transverse 
width increases much slower (like (logN) 13 , with (3 of the order of 1/3 for the particular 
values of parameters investigated numerically). We lack an analytical explanation of those 
findings. They show however that particles get more concentrated as N increases, suggesting 
that mean field treatment of their interactions should be relevant. 

In conclusion, we have shown how the physically motivated assumption of Bose condensa- 
tion of Lagrangian particles provides new predictions for the scaling properties of zero modes 
at large order and a clear-cut connection between zero modes and self-similar instantons at 
the same time. Not surprisingly, the Ctv curve shown in fig. [2] is not quite satisfactory on the 
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FIG. 2. Numerical estimates of scaling exponents (jv (square dots) versus IniV, obtained for N 
ranging between 1 and 10 6 , with 7 = 4/3 and d = 2. Data are almost perfectly reproduced by the 
analytical expression £ n = 2.87n/(1.0 + 0.7 , rT^pirT^iTi ) (solid line). 



low A/ side : in particular, it crosses the dimensional estimate (n = \- for iV ~ 10 rather 
than for N = 2 as it should [1|. It would therefore be useful to see how quadratic fluctuations 
(Bogoliubov approximation in the context of Bose liquid theory) correct the results obtained 
so far and possibly make them more accurate even at moderate values of N. Note that the 
inclusion of those fluctuations will restore some important symmetries in the ground state 
wave function, like rotational invariance. Preliminary steps in this direction suggest that 
this can be done rather harmlessly. Beyond these quantitative issues, an interesting outcome 
of this work is that it sets back on stage the instanton approach as a sound and potentially 
promising way for capturing, one day, the full complexity of Navier-Stokes turbulence. 
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